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Abstract 

We investigate the production of highly energetic top-quark pairs at hadron collid- 
ers, focusing on the case where the invariant mass of the pair is much larger than the 
mass of the top quark. In particular, we set up a factorization formalism appropriate 
for describing the differential partonic cross section in the double soft and small-mass 
limit, and explain how to resum simultaneously logarithmic corrections arising from soft 
gluon emission and from the ratio of the pair-invariant mass to that of the top quark to 
next-to-next-to-leading logarithmic accuracy. We explore the implications of our results 
on approximate next-to-next-to-leading order formulas for the differential cross section 
in the soft limit, pointing out that they offer a simplified calculational procedure for 
determining the currently unknown delta-function terms in the limit of high invariant 
mass. 



1 Introduction 



Top-quark pair production plays an important role in the physics programs of hadron colliders 
such as the Tevatron and the LHC. While much information about the top quark is already 
available, its high production rate at the LHC will eventually bring studies of its properties 
into the realm of precision physics. An especially useful observable is the differential cross 
section at large values of the top-pair invariant mass M. Many models of physics beyond the 
Standard Model predict the existence of new particles which decay into energetic top quarks 
and whose characteristic signal would be either resonant bumps or more subtle distortions in 
the high invariant mass region of the differential distribution. In fact, Tevatron measurements 
of the forward-backward asymmetry at high values of the pair invariant mass may already 
hint at the existence of such particles [UE]- Therefore, precision calculations of the pair 
invariant-mass distribution within the Standard Model are well motivated. 

The starting point for any study of two-particle inclusive differential cross sections such 
as the pair invariant-mass distribution is the next-to-leading order (NLO) calculations carried 
out roughly two decades ago [3J. Recently, these have been supplemented with soft-gluon 
resummation at next-to-next-to- leading logarithmic (NNLL) accuracy in [4,5j, building on the 
next-to-leading logarithmic (NLL) results of [6]. While such resummed calculations contain 
what are argued to be the dominant perturbative corrections at next-to-next-to-leading order 
(NNLO) and beyond, they suffer from two potential shortcomings. First, while they fully 
determine the coefficients of a subset of logarithmic plus-distribution corrections which become 
large in the soft limit z = M 2 /s — > 1, with the partonic center-of-mass energy, they 
do not fully determine the delta-function corrections in this limit at NNLO. The numerical 
contribution of these unknown coefficients as well as less singular terms in the soft limit are 
typically estimated through the method of scale variations, but this is by no means a fail-safe 
technique and additional information about the structure of these terms is valuable. Second, 
while [5] uses the parametric counting M ~ m t , when the top quarks are truly very boosted 
this counting breaks down. One must assume instead that M 3> m t , and recognize that 
resummed perturbation theory should also take into account powers of mass logarithms of the 
ratio m t /M. 

The primary aim of this work is to develop the theoretical framework needed to describe the 
invariant-mass distribution in the double soft and small- mass limit, and to explore some of its 
implications for perturbative predictions for the large-M distribution. The basic idea behind 
our approach is to weave together current understanding of factorization in either the small- 
mass or the soft limit into a unified description encompassing both. The first component is 
the factorization of partonic cross sections for highly boosted heavy-quark production worked 
out in [TJ. In the case at hand, where m t <C M, the results of that work imply that the 
partonic cross section can be factorized into a convolution of two functions: the cross section 
for massless quark production, and a convolution perturbative fragmentation functions for each 
of the heavy quarks. Given this form of the cross section for the small- mass limit, it is then an 
easy matter to perform the additional layer of factorization for the soft limit on the component 
parts. On the one hand, the massless partonic cross section in this limit can be factorized into 
a product of soft and hard functions using the techniques from [8], and on the other hand, the 
fragmentation function can be factorized into a product of collinear and soft-collinear functions 
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using the results of [9HT2]. A fully resummed cross section appropriate for both limits is then 
obtained by deriving and solving the renormalization-group (RG) evolution equations for the 
different functions separately. The anomalous dimensions appearing in the RG equations are 
known to the level sufficient for resummation of both mass and soft logarithms to NNLL 
accuracy. As simple and obvious as this approach is, it has yet to be fully worked out for any 
particular observable in top-quark pair production at hadron collider experiments. 

This formalism for the simultaneous resummation of soft and mass logarithms in the 
invariant-mass distribution is interesting in its own right. Moreover, with use of a proper 
matching procedure, it provides supplemental information to the current state-of-the-art pre- 
dictions based on soft-gluon resummation with the counting m t ~ M [5]. Particularly impor- 
tant in this regard is its use as a tool to calculate, up to easily quantifiable power corrections 
in m t /M, the full NNLO corrections to the massive hard and soft functions. Together, these 
pieces determine the coefficient of the delta-function coefficient in the fixed-order expansion 
at NNLO, a missing piece in currently available "approximate NNLO" formulas for generic 
values of the top-quark mass. Using our factorization formula for the double soft and small- 
mass limit, we can calculate the pieces of the NNLO delta- function correction enhanced by 
logarithms of the ratio mt/M. Furthermore, using the explicit NNLO results for the heavy- 
quark fragmentation function [13] and the virtual corrections to massless qq — > q'q' [H] and 
99 ~~ QQ [IS] scattering, we can very nearly determine the piece of the delta-function coeffi- 
cients which is constant in the limit m t /M — > 0. The missing piece is the NNLO soft function 
for massless partons, related to double real emission for gg — > qq and qq — > qq scattering 
in the soft limit. We do not calculate this function here, but plan to return to it in future 
work. While these delta-function pieces of the NNLO partonic cross section are of N 3 LL in 
the counting of soft-gluon resummation, including them can only make the predictions more 
precise and potentially strengthen the arguments in favor of the logarithmic counting underly- 
ing approximate NNLO formulas. This is currently an open point in soft-gluon resummation 
for the invariant-mass distribution, where the assumed dominance of logarithmic corrections 
is justified mainly through numerical studies at NLO and arguments based on dynamical 
threshold enhancement [16]. 

The remainder of the paper is organized as follows. First, in Section El we derive a fac- 
torization formula for the partonic cross section valid in the double small-mass and soft limit. 
This fixed-order expression contains large logarithms for any choice of the factorization scale. 
We deal with this in Section [3] by deriving and solving the RG equations for the compo- 
nent functions of the factorization formula, presenting in addition explicit perturbative results 
for their fixed-order expansions. In Section H] we combine those results into an expression 
for the resummed partonic cross section at NNLL in perturbation theory, and discuss differ- 
ent matching procedures needed to take into account power-suppressed terms away from the 
double small-mass and soft limit. In that section we also discuss approximate NNLO imple- 
mentations of the NNLL formula. Finally, in Section [5], we make preliminary explorations into 
phenomenological consequences of our results. We conclude in Section El 
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2 Factorization in the soft and small-mass limits 

We study the top-quark pair production process 

iVx(Pi) + N 2 (P 2 ) t(p 3 ) + t(p A ) + X(p x ) , (1) 

where Ni and N 2 are the colliding protons (or proton-antiproton pair), X is an inclusive 
hadronic state, and the top quarks are treated as on-shell particles. Two partonic channels 
contribute at lowest order in perturbation theory: the quark annihilation channel 

q(pi) + q(p2) -> t{ P3 ) + t{ P4 ) , (2) 

and the gluon fusion channel 

g(pi) + g{p2)-n(p3) + t(p4). (3) 

The momenta of the incoming partons are related to the hadron momenta by the relation 
Pi = XiPi (i = 1, 2). The relevant Mandelstam invariants are defined as 

s = (Pi + P 2 f , s = ( Pl +p 2 ) 2 , M 2 = (ps+prf , 

h = (pi - P3) 2 - ni 2 , m = (p2 - p 3 ) 2 - "i? • (4) 

In order to describe the invariant-mass distribution near the partonic threshold, it is convenient 
to introduce the following variables: 



M 2 M 2 n / 4m, 2 4m^ 

The quantity /3 4 is the 3- velocity of the top quarks in the ti rest frame. In the soft limit z — > 1, 
one has /5 — >- /5 t . Moreover, in that limit the scattering angle 9 is related to the Mandelstam 
variables according to 

M 2 M 2 
*! = -— (l-ftcosfl), ui = -— (l + ftcosfl), (6) 

and M 2 + t x + u\ = can be used to eliminate Wi as an independent variable. 

We will be interested in the double differential cross section with respect to the invariant 
mass of the top-quark pair and the scattering angle 9 in the parton center-of-mass frame. 
According to factorization in QCD, the double differential cross section can be written as a 
convolution of a partonic cross section with parton distribution functions (PDFs). We write 
this as 



d 2 a 8%(3 t 



f 1 dz 

^2 ~ ffij(T/z,iMf)Cij(z,M,nit,coa6,Hf) , (7) 



dMd cos 9 3sM 

where ///is the factorization scale. The parton luminosity functions functions Jfy are defined 
convolution of PDFs: 

ffijiv, v-f) = fi/m (y, fj>f) ® fj/N 2 (y, n/) ■ (8) 
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Here and throughout the paper the symbol <S> denotes the following convolution between two 
functions 



f 1 dr 

f(z)®g(z) = / — f(x)g(z/x). (9) 

J z X 

When there are several arguments in the functions / and g, the convolution is always over the 
first argument. As described in more detail below, we choose to define the PDFs with n\ = 5 
active light flavors, so that all physics associated with the scale of the top-quark is absorbed 
into the perturbative coefficient functions CV,-. These coefficient functions are proportional 
to differential partonic cross sections. Our aim is to study the factorization properties of 
these partonic cross sections in the double soft and small- mass limit, where (1 - z) < 1 and 
Our strategy is to first discuss the soft and small-mass limits separately, and then 
combine them into a single formula which is true for both limits simultaneously. 

Factorization of differential partonic cross sections in the soft limit has been studied in 
[5|l8l fT71 - [25] . In the soft limit z — > 1, the partonic cross section can be factorized into a hard 
function and a soft function according to 

djiz, M, m t , cos 9,fi f ) = Tr [h%(M, m t , t 1}f ji f ) - z), m t , t uf i f )j + 0(1 - z) , (10) 

where we have used that in the soft limit dependence on the scattering angle can be expressed 
in terms of tx, see The superscript m on the hard function HJ? and the soft function 
S?j indicates that they are computed with finite top-quark mass, as opposed to the massless 
hard and soft functions introduced below. Both of these functions are matrices in the space 
of color-singlet operators for Born-level production. The hard function is related to virtual 
corrections to the two-to-two scattering processes qq(gg) — > tt. The soft function is related to 
real emissions in the soft limit, or more precisely to the vacuum expectation value of a Wilson- 
loop operator built from time and light-like Wilson lines associated with soft emissions from 
the heavy and light quarks. Note that in this limit, we only need to consider ij = qq, gg, since 
the qg channel is suppressed by powers of (1 — z). 

Factorization of partonic cross sections for heavy-quark production in the small-mass limit 
was considered in [7j . It was shown that partonic cross sections in this limit can be factorized 
into a product of massless cross sections with perturbatively calculable heavy-quark fragmen- 
tation functions. Generically, for a cross section differential in the energy fraction z = E/E max 
of the top quark, this factorization is written as (see, for instance, [T3] ) 



a J z 

where da a /dx is the MS-renormalized differential cross section for the production of a massless 
parton a, and D^fi is the heavy-quark fragmentation function defined using an a s with nj 
active flavors. The sum over the massless partons labeled by a includes the case a = t, and 
the heavy quark is considered massless in the calculation of da t . Much as the PDFs describe 



^^More precisely, we work in the limit where the Mandelstam variables s, ti,ui 3> m t . 
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radiation collinear to initial-state partons, the heavy-quark fragmentation functions describe 
radiation collinear to the energetic final-state heavy quarks. The heavy-quark fragmentation 
functions are however perturbatively calculable, since the mass of the top quark serves as a 
collinear regulator. 

In the case of the invariant-mass distribution in top-quark pair production at hadron col- 
liders, we need to modify the generic formula ( II ip in several ways. First, since this observable 
contains information about both the top and the anti-top quark, we need a fragmentation func- 
tion for each of them. Second, we must introduce heavy-flavor coefficients related to matching 
six-flavor PDFs onto five-flavor ones, which induces an additional source of mt-dependence 
into the formula. Finally, although not strictly necessary, we will follow [12] and also perform 
such a matching for the fragmentation functions. The matching relations between the PDFs 
and fragmentation functions in the n\ + n^- and nj-flavor theories are 

D% l t +nh \z,m t ,Lif) = C a/t (z,m t ,Lif) ® D^\z,m t , fi f ) , (12) 
ff% l+nh \z,m hP L f ) =C${z,m t ,vi f )®jF$ i \z,iM f ). (13) 

The heavy-flavor matching coefficients C a /t and C 1 ^ on the right-hand side of the above equa- 
tion are proportional to powers of rih — 1. They are obtained by comparing partonic matrix 
elements with and without the top quark as an active flavor, and are known to NNLO in fixed- 
order for both the fragmentation functions [T2] and the parton luminosity functions We 
will encounter them again in Section HJ when we discuss the RG running of massless coefficient 
functions to scales below the flavor threshold at fit ~ m t . 

Taking these points into account, the factorization formula for the coefficient function (J2J) 
in the small-mass limit readsEI 

dj (z, M, mt, cos 6, pL f ) = C#(z> M, h, ///) <8> B^\z, m u N ) 

a.b 

® C a/t (z, m t , [if) <g> C b/ t(z, m t ,ii f ) ® C ff (z, m u ji f ) + O > ( 14 ) 

where the sum is over all parton species a,b G {t, i, q, q, g}. The functions are the partonic 

cross sections obtained from the massless scattering process ij — » ab + X, where X indicates 
additional final-state partons. The objects are defined as the following convolution of 

heavy- quark fragmentation functions: 

B^\z,m t ,fi f ) = D^(z,m t ,fi f ) ® D^\z,m u fi f ) . (15) 

This convolution of heavy-quark fragmentation functions is completely analogous to that defin- 
ing the parton luminosities in (jSJ). It arises after generalizing (1111) to a two-fold convolution 
and performing a change of variables. 

2 We have used a slight abuse of notation and expressed the dependence on the scattering angle in terms of 
t\, which is in general only possible in the soft limit. In converting to the scattering angle, we keep the exact 
mass dependence in t\ according to otherwise a t-channel singularity emerges upon integration over 9. 
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We are now ready to discuss the joint limit z — > 1 and m t /M — > 0, which is the main 
theme of this paper. The key point is that these two limits are independent and commutative, 
so that we can take them one- by-one in any order and obtain the same result. We choose 
to start from the factorization formula (Tl4|) for the small- mass limit, and then study the 
behavior of its component parts in the limit z — > 1. We thus discuss the factorization of the 
massless coefficient functions and the fragmentation function in the soft limit. The alternate 
method of starting from the factorization formula ffTOj) for the soft limit and then studying 
the factorization of its component parts in the small-mass limit is discussed in Appendix |A] 

We first deal with the massless coefficient function Cf^. To factorize it in the soft limit, 
we observe that nothing in the derivation of factorization for the massive coefficient function 
( TTUj) makes reference to the mass of the top-quark. Therefore, the form of factorization for 
the massless coefficient function is exactly the same. The result is thus 

CffaM^Hf) =Tr [iZ^M, -z)M,Hf)] +0(1-*). (16) 

We have used that only a = t contributes to (TT2"|) at leading power in (1 — z). The hard 
function Hij is obtained from virtual corrections to two-to-two scattering with massless top 
quarks, and the soft function Sij involves only light-like Wilson lines related to real emission 
from massless partons. The top quark is treated as massless in both the external states and 
in internal fermion loops, so both the hard and soft function are defined in a theory with six 
active massless flavors. 

The factorization of the fragmentation functions in the z — > 1 limit was explained in p-TTT]. 
and also within an effective field-theory framework in [12]. The main result of those works 
is that after the matching onto the n^-flavor theory as in (fl2j) . the fragmentation function 
factorizes into a product of two functions: one depending on the collinear scale m t , and the 
other on the soft-collinear scale m t {\ — z). We write this factorization as 

D$(z,mt,Hf) = C D (mt,fi f ) S D (m t (l - z),fXf) + 0(1 - z) . (17) 

The fragmentation of t to t follows the same factorization with the same coefficient functions. 
The soft function So is related to soft-collinear emission and is equivalent to the partonic 
shape-function appearing in 5-meson decays [TTjlT2"]. The matching coefficient Cd is indepen- 
dent of z and is a simple function related to virtual corrections. 

Combining all of the information above, the factorization formula for the partonic cross 
sections in the joint soft and small-mass limit is 

Cij{z,M,mt,coad,Hf) = C 2 D {m u pf) Tr Hij(M, t u nf) S i:j {\fs{l - z),t u Vf) 

® Cjj(z, m t ,fi f )0 C t/t (z, m u n f )® C t/t (z, m t ,fi f ) 

® S D (m t (l -z),n f )® S D (m t (l - z), n f ) + 0(1 - z) + O (^) . (18) 

The factorization formula flTSJ) is the central result of this section. In the limit in which it is 
derived, any choice of \ij generates large logarithms in the soft or small-mass limits. We deal 
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with this problem in the next section using RG techniques. In deriving and solving the RG 
equations it will be useful to introduce the Laplace-transformed functions 



Cij(N, M, m t , cos 9, /if) 



M 2 

Sij \ ivv' ''^ 



c t J I ^j^,m t ,(i f 



d^e ^ N Cij (z, M, m t , cos 9, /i/) 

POO 

Jo 



d^e iN &l(z,m u Hf) <g) C t/t (z,m t , /if) ® C t/t (z,m t , jif] 



(19) 



where £ = (1 — z)/y/z and N = Ne lE . In Laplace space, the factorization formula becomes a 
simple product of the different functions and reads 



c ij (N,M,m t ,cos9,fi f ) = C 2 D im u fi f ) Tr 



Hij(M,ti,fj, f ) Sij In 



M 2 



cm± 



(20) 



3 The matching functions: fixed-order expansions and 
RG evolution 

The component parts of the factorization formula (fl8l) can be viewed as matching functions 
in effective theory. In this section we explain the one-scale calculations needed to extract 
the fixed-order expansions of these matching functions, present their RG equations and the 
solutions thereof, and give the ingredients needed to evaluate these RG-improved matching 
coefficients at NNLL. While such NNLL calculations require only the NLO perturbative ex- 
pansion of the matching coefficients, we also collect all current knowledge at NNLO, part of 
which we will use in our numerical analysis later on. 



3.1 Hard function 

The hard function is related to virtual corrections to the two-to-two scattering processes 
underlying Born-level production. The method for calculating the hard-function matrix for 
the counting M ~ m t in fixed-order perturbation theory was described in detail in [5], where 
results valid to NLO were given. This boiled down to calculating color-decomposed UV- 
renormalized on-shell scattering amplitudes and subtracting poles in the 4 — d = 2e — > limit 
using an IR renormalization factor. 
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There are in fact two ways to calculate the hard function in the massless case. The first is to 
set mt = at the start of the calculation and follow the same procedure as for the massive case. 
Then the UV-renormalized scattering amplitudes and the IR renormalization factors change 
compared to the massive case, but the method for extracting the hard function is exactly the 
same. The second way is to Taylor-expand the massive result [5] in the limit mt — > and 
then convert the result to the massless case using the relation between massive and massless 
amplitudes in the small-mass limit derived in [27] . We have calculated the NLO hard function 
using both methods and confirmed that they agree. This NLO result for massless scattering 
is not new: it is actually a special case of the more general results for four-parton scattering 
given in [28], using the one-loop calculations of [29]. We have checked that we reproduce those 
results using the procedure described below. 

As opposed to the massive case, where only a limited set of NNLO virtual corrections have 
been calculated [301437] , both the one- loop times one-loop and two- loop times Born interference 
terms are known for the case of massless two-to-two scattering [T ^[38TH0] . These results are 
implicitly summed over colors and cannot be used to extract the hard matrix directly, on 
the other hand they can be used to extract the contribution of the NNLO hard function 
to the approximate NNLO formulas covered below. Although we do not go through this 
straightforward but tedious exercise here, it is an important point that all of the diagrammatic 
calculations are in place. 

In order to give explicit results valid to NLO we define expansion coefficients of the hard 
function H as 



// 



8d 



R 



H (0) 



Air' 



Cts 

Air 



(21) 



where d^ = N in the quark annihilation channel and d^ = N 2 — 1 in the gluon fusion channel, 
with N = 3 colors in QCD. The LO result H® is trivially obtained from the formulas in [5], 
which is regular in the limit mt — > 0. For the qq channel we have 



H, 



(o) 



QQ 



0\ t\ 



2 



M 2 



(22) 



and the result for the gg channel is 



H (o) 

99 



( 1 
' N 2 



1 ti —Ul 



N 



\ 



M 2 
N 



1 tj— Ul 



N 

Ci 



M 2 
-«i) 2 



M 4 
ti—ui 
M 2 



ti—ui 
M 2 

1 J 



tj + U 2 

2tiUi 



(23) 



We do not list the explicit result for the NLO hard function H^ l \ since it has essentially been 
given in [28J. For what concerns the quark annihilation channel, we convert those results to 
our case by extracting the elements of the hard matrix from Eq. (39) of [28] . Subsequently, 
one needs to consider crossing symmetry and to permute the arguments of the various matrix 



elements according to Hij(s,t,u) — > 



Hij(u 



s,t), as explained in Table 1 of 



Finally, 



one should exchange the element indices 1 <H- 2 to match the notation employed here, and 
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compensate for an overall factor, so that 



1 



H 22 (u,s,t) H 21 (u,s,t) 
H 12 (u,s,t) H n (u,s,t) 



(24) 



where all of the elements of the matrix on the right-hand side are taken from Eq. (39) of [28] . 
The gluon fusion case is slightly more complicated because the authors of [28] use a different 
color basis. In that case, the NLO correction can obtained from Eqs. (56) in [28] through the 
rotation 



99 



AO 1 



NLO 



ttNLO 

NLO 
13 



ttNLO 
n 12 

ttNLO 
n 22 

ttNLO 
n 23 



ttNLO\ 

ttNLO 
n 23 

Hg L ° 



where the matrix elements on the right-hand side are from 



O 



and 



(25) 



O 



( J- I I\ 

' 2JV 2 2 4 

1 _ 1 1 

2JV 2 2 

Vi o o/ 



(26) 



In the resummed formulas below we will need an expression for the hard function evaluated 
at an arbitrary scale /if, given its value at an initial scale fi^ ~ M where it contains no large 
logarithms and the fixed-order expansions above can be applied. This is obtained by deriving 
and solving the RG equation. As with the matching coefficient itself, we can derive the RG 
equation either by setting m t = from the start and using the exact same methods as the 
massive case [5], or we can start from the massive result and take the limit m t — > using 
the relations between massive and massless amplitudes mentioned above. We have checked 
that the two methods agree. In any case, the RG equation for the massless hard function is 
completely analogous to the massive case studied in [5] and is given by (here and below we 
suppress dependence on the channels qq and gg when there is no potential for confusion) 

-^—H{M } tt, n) = T H (M, h, /x) H{M, t u At) + H{M, t u n) T j H (M, h,^. (27) 

The explicit result for the anomalous dimension matrix to two-loop order in the color basis 
of [5] is easily derived by making use of the general result [H]|l2] f° r massless scattering 
amplitudes, and reads 



2C F 7 cusp (a s ) In 



cusp 



at. 




m + 47 9 (a s ) 



+ 7cusp(« s ) In 4 2 \ 



{2i 



and 



99 



M 2 

(N + C F ) 7 cusp(« s ) ( In — - m ) + 2 1 9 {a s ) + 2 7 > s ) 
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+ iV7, 



cusp j 



111 



M 2 



in 



/o o o\ 

1 

\o o l) 



/o 



7cusp(*^s) 1^ o 



o \ 



1 TV iV 2 -4 

4 4/V 

V° ? -f/ 



(29) 



For convenience, we have collected results for the various scalar anomalous dimension functions 
in Appendix O 

Since the evolution equation has the same structure as in the massive case, it can be 
solved using the same methods. In presenting the solution it is convenient to decompose 
the anomalous dimension into a logarithmic piece multiplying the unit matrix and a non- 
logarithmic part containing the non-trivial matrix structure: 



T H (M,t 1 ,fi) = A{a s ) In 



M 2 
~7 



Z7T l+-y h (M,t u a s 



(30) 



where A = 2C F -f cusp = 2T q cusp in the qq channel and A = (N + C F )7 cusp = r^ usp + T q cusp in the 
gg channel. The solution to the RG equation can then be written as 

H(M, h, //) = U(M, ti,n h , n) H(M, ti,fM h ) C/ f (M, t U fi h , //) , 

with 



(31) 



U(M,ti,ii h ,ii) = exp 



M 2 

2S A (Vh, y) - a A (n h , y) ( In — - m 



u(M,ti,n h ,n) 



(32) 



The RG exponents are given by 
<Sa(/^, A 1 ) = - I da 



A{a) f da' 
W) J W) 



«s(/i) 



da 



A(a) 



(33) 



where (3(a s ) = da s /d\nfi is the QCD /^-function (whose expansion coefficients are given in 
Appendix [C]). The matrix- valued contribution to the evolution function reads 



u(M, h,fih, n) = Vexp 



da 



~f h (M, ti, a) 



(34) 



The exact solution to the RG equation is evaluated as a series in RG-improved perturbation 
theory as in [5], where explicit expressions valid to NNLL order were presented. Such an 
NNLL calculation requires the NLO corrections to both the hard matching function (a one- 
loop calculation) and the anomalous dimension (a two-loop calculation). 



3.2 Soft function 

The soft function is related to real emission corrections to massless qq, gg — > ti scattering in the 
soft limit. A more formal definition in terms of the vacuum expectation value of a Wilson-loop 
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operator was given for the massive case in [5]. The position-space result for this object can be 
directly converted into the Laplace-transformed function (TT9|) . We can adapt that definition 
to the massless case simply by changing time-like Wilson lines representing emissions from 
massive particles to light-like Wilson lines representing emissions from massless ones. We will 
present results for the Laplace-transformed soft function directly. We define the perturbative 
expansion of this function as 



;(0) 



4lX 



(i) 



Air) 



(35) 



At lowest order, the result depends on whether the initial-state partons are quarks or gluons, 
but not on the mass. The result is 



g (o) 




(36) 



in the quark annihilation channel, and 
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(37) 



in the gluon fusion channel. 

To obtain the Laplace-transformed soft function at NLO we evaluate the following position- 
space integrals [5] 



Tij{e,x Q ,n) = - 



(4:7111 



2\e 



7i 



2-e V i' V J 



d d k 



ik°xo 
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(2tt) 5(k 2 ) 9(k°) 



(31 



where are the light-like four- velocities of the partons from the Born-level scattering process. 
When i = j these integrals vanish since vf = 0, while for i ^ j they are equal to 
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(39) 



where 



In 



4 



(40) 



We make use of this result by expressing the scalar products V\ ■ Vj in terms of the Mandel- 
stam variables, by subtracting the IR poles in MS, and by Laplace-transforming the integrals 
through the replacement L — > —L. This leads to 



1 



12 



IT 



L z + — 
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Xi 3 = X 



21 



(L + ln(r)) 2 -^-2Li 2 (l-r) 
o 



Xu = X 23 = - (L + ln(l - r) f - ^ - 2 Li 2 (r) 

o 



X n = X : 



22 



X33 — X 



11 







(41) 



where r = —t\jM 2 . One obtains the matrix- valued soft function in Laplace space by evaluating 



(42) 



where the ly^ are the color matrices from [S], which are different for the gg and channels, 
but make no reference to the parton mass. 

To obtain the NNLO correction to the soft function requires a new calculation which is 
beyond the scope of this paper. However, as a compromise, we can use the RG equation below 
to derive all of the coefficients proportional to powers of logarithms in the Laplace-transformed 
NNLO correction, whose form is 



s^(L,M,t 1 ) = J2 s n(M,t 1 )L n 



(43) 



n=0 



The results for the coefficients are fairly lengthy and since we use them in this paper only for 
the factorization check described in Section l4~2l we do not list them here. 

The Laplace-transformed functions are the central objects used in solving the RG equations 
below. One can also convert them to the momentum-space functions using a set of replacement 
rules. In the case where the first argument of the soft function is expressed in terms of 
Vs(l — z) — M(l — z)/yfz, the resulting distributions are 



PL(z) 



M 2 (l-z) 
jj?z 



(44) 



As shown in [5] , the momentum-space soft function is derived from the Laplace-space function 
by making the replacements 

l->8(\-z), 

fM 2 \ 

I? ^ 4P[(z) + 8(1 -z) In 2 (*£Y 
L 3 ->■ QP^(z) - ^ 2 Pq(z) + 5(1 - z) 



ln d 



M' 
~J 2 



+ 4C 3 
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L 4 -> 8Pg(z) - 16n 2 P[(z) + 128( 3 Pq(z) + 5(1 - z) 



In 4 — + 16C 3 ln — 



In order to translate the PL into the conventional P„ distributions, 



M 



2\ 1 



(45) 



PJz) 



1 - z 



\n n (1 - z) 



(46) 



we employ the general relation 



fc=0 

fe-i 

+ £ 

j"=o 



2 J (-ir 



fe-i 



2 fc P fc (^) 



ln^l-z) \n k ~ j 
l-z 



(47) 



-5(l-z) 



1 , \n j (l-x)ln k ~ j x 

ax 
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The numerical affects of keeping the power-suppressed terms proportional to ln m z/(l — z) was 
discussed in detail in [51125]. 

The soft function obeys a non-local RG equation which is solved using the Laplace trans- 
form technique [13]. The RG invariance of the total partonic cross section implies that the 
Laplace-transformed soft function satisfies 
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<iln/i 



M 2 

s I In— , M,ti,/x 



M 



A(a s )ln— + 7 st (M,t 1; a s ) 



— s ( In M, t 1; /i 



M 2 

In— ,M,*i,/i 



M 2 



M 2 

A(a s )ln— + 7 s (M,t 1; a s ) 
/z 2 



We have defined 



Y(M,t u a s ) =~ f h (M,t 1 ,a s )+ [2 7 *(a s ) + 2 7 *«(a s )] 1. 



(48) 



(49) 



Note that while the form of the anomalous dimension is analogous to that in the massive 
case, it picks up an extra term 2^ q , which is needed to cancel the /x- dependence from the 
fragmentation function as determined by the RG equation (I57)) below. 
The solution for the momentum-space soft function reads 



5(o;,M,ti,/i / ) = a/I exp [-AS A (/i s , fi f ) + 4ay(/x s , /i/) + 4ay, (/x s , ^ 



X W T (M,ti,/i f ,/i s ) s(dr, A ,M,ti,ljL s ) u(M,t l7 fif,fi a ) - — 



CO 



r(277 A ) 



(50) 



where one is to set t/a = 2a A (ii s , /if) after performing the derivatives. For values of 2rjA < 0, 
the w-dependence must be interpreted in the sense of distributions. 
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As always in RG-improved perturbation theory, the aim of a formula such as (loPl) is to 
allow one to evaluate the soft function at an arbitrary scale fif given its result at a scale /x s 
where it is free of large logarithms. However, the question of what exactly this fi s should be 
is currently a source of debate in the literature. We will discuss this issue in more detail when 
presenting results for the RG-improved partonic cross section in Section HI 



3.3 Fragmentation function 

The perturbative fragmentation function was calculated at NNLO for generic values of z in 
In this section we focus on the parts of that result required for the resummed analysis, namely 
the leading terms in the soft limit of the function with n\ active flavors defined in ffT2l) . In 
particular, we list results for the functions Sd and Cd appearing in the factorized form (TTTT) . 
determined previously in [T2j . 

The function Sd is related to the partonic shape-function in 5-meson decays and can be 
derived from the two-loop calculations in [H] . We define its perturbative expansion in Laplace 
space as 

8D = 1+ 4^* D + D + '"' (51) 

The expansion coefficients with N = 3 colors are 

IOtt 2 
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(52) 



104 76vr 2 16 
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— rii \ L A + 1 1 rii L 

27 7 V 9 27 27 

16 16tt 2 

" 72Ca + 81^ " ~rf ni 

12137T 4 1132C 3 , / 16 147T 2 
H T77z h -— + + nj. (53) 



405 9 V 243 27 27 

It is a non-trivial check on the factorization formula (|17[) that the plus-distributions in the 
fragmentation function are all related to the momentum-space representation of this function, 
obtained by the set of replacement rules analogous to (145|) . but with the substitution M(l — 
z)/yfz -> m t (l - z). 

The coefficient Cd is related to virtual corrections to the fragmentation function and is a 
simple function independent of z. Since the NNLO correction to the fragmentation function 
obtained in [TJ] was not split into real and virtual corrections, it is not possible to obtain the 
coefficient directly from that work. Instead, it must be determined by using the result for the 
shape-function given above along with the factorization formula ()17|) for the fragmentation 
function. Defining expansion coefficients in analogy to fl5~T|) . we have (with N = 3 colors) 

C$> (m t , fi) = i (L 2 m + L m + 4 + ^ V (54) 
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(2) x_8 r4 (20 8 \ 3 /406 28tH 52 \ 2 



^'(m, /•) = 5 C + ^ - ^ j Li + ^- - — - -n, j L 
/2594 248tt 2 232C 3 308 16tt 2 \ 

+ {-w + —7 3--^r ni -^r ni ) Lm 

21553 107vr 2 749vr 4 260C 3 16vr 2 , /1541 74tt 2 104C 3 \ 

H 1 1 1 In 2 — 1 1 m , 

162 3 405 9 9 \ 243 81 27 J 

(55) 

with L m = ln(/i 2 /m 2 ). We discuss a possible cross-check of this result in Appendix fAI 

In the factorization formula for the invariant mass distribution we need the convolution of 
two fragmentation functions, which up to NNLO has the form 

M(z, m t , n) = 8{1 - z) + 2 (g) D^(z, m t , /x) 

+ (g) 2 [2D< 2 )( Z) m 1 , /1 ) + fl( 1 )( Z , ffll ,^D( 1 >(z, mi ,ri] . (56) 



In order to evaluate the convolutions between the different plus distributions in the last term 
we use the formulas derived in Appendix iBl 

The RG equation for the fragmentation function is a non-local one. It is given by 



( t/ £(z,m t ,n) = P qq (z,ii) ® D\i%{z,mt,n) ( 57 ) 



where P 9g is a time-like Altarelli-Parisi splitting function whose structural form in the soft 
limit is 

2T q (a ) 

P gq (z, pi) = 7 s) + 2rfr{a.)8{\ - z) . (58) 
U _ z )+ 

From this equation, and the fact the So is equivalent to the perturbative shape-function from 
5-meson decays [TTl[T2"] . the RG equations for the function Cd can be derived. The RG 
equation for Co is local, while that for Sp is non-local and solved using the Laplace transform 
technique. The result for the RG-improved fragmentation function reads 



D(z, m t ,fi f ) = exp [2S r ^ cuBp (fi ds , fi dh ) + 2a 7 s(/i ds , fi dh ) + 2oy, (fi f , fi dh )] ( — j 



-~/EVd / m \ T >d 



-2a T q (jif,fldh) 



C D {m^ dh )~s D {d m ^ ds )^^—j (1 _ ^ , (59) 

with r] d = 2a r 9 usp (/ij, pi ds ). The explicit results for the anomalous dimension 7 5 can be found 
in [12] as well as in Appendix [C] of this paper. While it is obvious that the scale choice 
(J'dh ~ m t eliminates large logarithms in the coefficient function Cd, the choice of the scale fi ds 
is again a debatable point which we come back to later on. 
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3.4 The heavy-flavor coefficients 



The definition of the heavy-flavor coefficients C t /t and C 1 ^ was given in (|T2|) . The partonic 
matrix elements needed to evaluate those expressions to NNLO are known from [13] for the 
fragmentation functions and for the PDFs. Rather than give the results separately, we 
quote only the result for the Laplace-transformed combination of the three functions appearing 
T9l) . Expressed in terms of a s with five active flavors, the result is 
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(60) 
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(61) 



with 



If desired, these can be converted to momentum space by the set of replacements 
H = M. 

The RG equations for the heavy-flavor coefficients Cff and Ct/t follow from the fact that the 
parton luminosity and fragmentation functions in the rih + n\ flavor theory obey the standard 
Altarelli-Parisi equations in the soft limit. This implies 



d 



(iln/x 



C t/t {z, m t , aO = [P^ +ni {z, fi) - P£(z, fi)] (8) C t/t (z, m t , fj) , 



(62) 



(iln/i 



Cj f (z,m t ,fi) 



2 [P^ +ni {z^) -Pg(z,n)] ®C}(z,m t ,ri 



(63) 



where we have used that non-diagonal evolution is suppressed by powers of (1 — z). The 
function P gg is defined in analogy with (158]) after the obvious replacements, and both T^ usp = 
CaIcusp and 7^ 9 can be read off to two loops from Appendix O The superscripts indicate the 
number of active flavors to be used in both a s and the coefficients of the anomalous dimen- 
sions themselves. When expressed in terms of a common five-flavor coupling, the anomalous 
dimensions P™ h+ni pick up explicit powers of lnrrit/n related to the a s decoupling relation 



a 



{n h +ni) 



a 



(ni) 



2 r at ] 
+ nh 3 m ^T + 



(64) 



which implies that the anomalous dimension for the heavy-flavor coefficient no longer has the 
simple form ( |58|) . We note that for the NNLL analysis, we need only the NLO correction from 
cf 9 , which contains no large logarithms for fi t ~ mt- Beyond NNLL accuracy, however, such 
a choice leads to large logarithms in 1 — z that are not resummed. On the other hand, since 
the scale for logarithms related to the decoupling relation fl6"4"l) is unambiguously ji ~ m t , it 



16 



is conceivable that one could derive a method for resumming logarithms between the scales 
m t (l — z) and m t within the heavy-flavor coefficients. We will however leave this as an open 
point in our analysis of resummation in the double soft and small-mass limit. 



4 The partonic cross section at NNLL and approximate 
NNLO 

In this section we derive the final expression for the resummed partonic cross section at NNLL 
in the double soft and small-mass limit. We discuss a few points having to do with its practical 
implementation, and then turn to its approximate NNLO implementation. 

4.1 Partonic cross section at NNLL 

To derive the final result for the resummed partonic cross section at a scale fif we insert the 
RG-improved results for the hard, soft, and fragmentation functions presented above into the 
factorization formula (II 8p . The convolution integrals can be performed analytically using 

1 dz > 1 I „ F fa) r ^) 1 (65 ) 

(i - z') 1 ^ (i - z/z') 1 -^ rfa + r)2) (i - zy-m-^ ' 

where the approximation is true in the limit z — > 1. We can also simplify the various products 
of evolution matrices by employing the relations 

u (M, cos 9, fif, Us) u (M, cos 9, fih, /-*/ ) = u (M, cos 9, fih, H> s ) > 

dAiPs, fJ>h) + <U(A%, At/) = CL A ((A 8 , (if) , 

Sa (Hh,Hf) - S A {y>s, Hf) = S A {Hh, V>s) - a A {n s ,Hf) In — . (66) 

The only subtlety is related to the treatment of heavy-quark threshold effects at fi t ~ fn t . 
However, this is a standard problem in RG-improved perturbation theory involving heavy 
quarks, and we deal with it in the usual way [15]. To understand the logic, it suffices to consider 
a hypothetical observable involving three widely separated scales M ^> m t 3> /Uo, which 
satisfies a factorization formula of the form C(M, fi)D(m t , fi)F(fi Q , fi). If C is a coefficient 
function whose RG running is known, and the goal is to evolve it from a high scale /im ~ M 
to the scale /xo below the heavy-flavor threshold, one uses the following schematic equation 

C^\M,m t ,fi ) = U^\M,^ii t )M h (m t ,ii t )U^ +n *\M,^ . (67) 

In words, one uses six-flavor evolution functions U (jll+nt ^ above the flavor threshold at /i t , and 
five-flavor evolution functions below it. The change in the number of flavors induces a 
matching coefficient at the flavor threshold. It is determined by requiring that the partonic 
cross section in the rif and rif — 1 flavor theories be equal at the scale [it. 

C {n f\M,fi t )(DF)^ = & n t- 1 \M,7n t ,lh){DF)( n t- 1 ') , (68) 
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where the ()w denotes a partonic matrix element in the theory with / massless flavors eval- 
uated at the scale fi t . 

The generalization of this simple picture to our case is straightforward. The role of the 
coefficient C is played by the massless hard and soft functions, and that of the heavy-flavor 
matching coefficient by the convolution of the three functions appearing in the second line 
of (fT8"j) . The explicit result reads 



C(z, M, m t , cos 9, fif) = exp 



4^r« usp (/Ms, Hdh) + 4oy>(/i t , fif) + 4cy, 9 (fi t , Hdh) + 4a T s(/i fis , fi dh ) 

exp[4a 7 «i>(/i s , fi t ) + 4a 7 0, (/x s , //t)] n/ =e 
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+ 2a r?usp (^,/^ s )ln-^ 



n.f =5 



x Tr 



X s 



U{M, t u iM h , fi s ) H(M, h, fi h ) U\M, t u fi h , fi s ) 
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ln — + d n >,M,t 1 ,n, 

Ms 



Ct(d v >,m t ,iJt) 



In-r- + d v ',m t ,fi dh ,IJ.ds 



-2<y B v' 



f=5 r(2 V >) (i - ^)i-v 



(69) 



and in addition 



We have defined 77' = [2a j4 (/x s , ^)] n/=6 + [2a A (fi t , fi f ) + 2a r?u8p (/i / , /i ds )] 

d(L,m t , fidh, Hds) = [Co(^t) ^dh)Sn(L/2, fids)} 2 - We have indicated with the subscripts the 
number of active massless flavors rif to be used in evaluating the running coupling constant 
and perturbative functions in the various parts of the formula. This number of active flavors 
is chosen according to the physical picture of the schematic example (j67j) above, namely that 
of integrating out heavy degrees of freedom until reaching a scale under the flavor threshold 
below which the remaining degrees of freedom are factorized into the PDFs. However, it is 
formally true for any value of the factorization scale. It thus provides a convenient way to use 
the standard PDFs with five light flavors even when ///is far above the heavy-flavor threshold 
at fit ~ wit, as it explicitly resums any large logs in mt/M in the partonic cross section for 
such a scale choice. 

The result (169]) is the final expression for the resummed partonic cross section in momentum 
space. It can be evaluated perturbatively at NNLL order using the results given in the previous 
section. There are two important issues in terms of its numerical evaluation. The first is a 
technical one having to do with the choice of matching scales, the second is a practical one 
having to do with the power corrections away from the soft and small-mass limits. We end 
this section on the resummed cross section by discussing these in turn. 

As alluded to several times in the previous section, the philosophy of RG-improved pertur- 
bation theory is to use RG evolution factors to evaluate the matching functions at an arbitrary 
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scale /jbj given their value at an initial scale where they do not involve large logarithms. This 
RG running then exponentiates large corrections appearing when ///is parametrically far from 
the natural scale. For the hard function and the coefficient Cd, which obey local RG equa- 
tions, the choice is straightforward: one uses /i^ ~ M and /t^ ~ m t- For the massless soft 
function and the function So, the correct choice of this scale is less obvious. If the goal is to 
resum logarithms of (1 — z) in the partonic cross sections, then the natural scales are M(l — z) 
and m,t(l — z). However, such choices are ill-defined at the level of the momentum-space 
result fl69|) . Partonic logarithms can be resummed by choosing the scales /i s and fids at the 
level of the Laplace-transformed functions ( 119]) and performing the inverse transform back 
to momentum space numerically, but at the cost of introducing the Landau-pole singularity 
familiar from Mellin-space implementations of soft-gluon resummation for top-quark pair pro- 
duction [1&H52]. An alternate method is to choose the two scales as numerical functions of M, 
in such a way that the logarithmic corrections to the hadronic cross section arising from those 
in the partonic one are minimized after convolutions with the PDFs [16]. The fixed-order 
expansions of the resulting expressions at any finite order in the logarithmic counting are then 
of a different structure than those in the partonic cross section [2o"|lo*3] . Studying the numerical 
differences between these methods would be an interesting exercise, but since we will not do 
detailed phenomenology in the current paper we leave this issue aside. 

Dealing with the power corrections away from the double soft and small-mass limit is also 
important, although considerably more straightforward technically. The standard method 
is to include these corrections at NLO in fixed-order perturbation theory, thus obtaining 
NLO+NNLL accuracy. This is accomplished by evaluating partonic cross sections da as 



da(n f ) = da NNLL ({/ii}) 



+ ( da NW ^ f ) - d<7 ({/tj} = ///}) 



(70) 



mt— >0,z— >1 



where by {/i,;} we mean the set of scales /t?j,/t s , ... appearing in (1691 . Such a resumma- 
tion formula is useful only at values of the invariant mass where m t /M is truly small. It is 
straightforward to extend the matching procedure to take into account in addition the set of 
higher-order m t /M corrections determined from soft gluon resummation at NNLL with the 
counting m t ~ M used in [5], thus yielding a result useful for the full range of M, but we do 
not give the explicit results here. 



4.2 Partonic cross section at approximate NNLO 

The resummed formula from the previous subsection can be used as a means of extracting what 
can be argued to be dominant part of the full NNLO correction in fixed-order perturbation 
theory. This truncation of the resummed expansion is of course not valid if the logarithms are 
truly large, but it is easy to imagine a situation where (1 — z) and m t /M are good expansion 
parameters for the fixed-order corrections, but not necessary so small that logarithmic correc- 
tions beyond NNLO are numerically important. In this section we focus on such approximate 
NNLO formulas based on our NNLL results, and explain how the results in this paper can 
offer an improvement on those previously derived in [HE]. 
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It is convenient to discuss the approximate NNLO formulas at the level of Laplace-transformed 
coefficients. The general expression for the NNLO correction in Laplace space reads 



c(N, M, m t , cos 9, jif) 



g (0) (M,m t , cos 9, fi f ) + (^j c {1 \N,M,m t ,cos9,n f ) 



2 c (2) (N, M, m t , cos 9, fi f ) + 0{a\ 

4:71/ 



c {2) (N,M,m t , cos 9, fi f ) = 



n=0 



M 2 ( 1 

^\M,m t ,t x , N ) ^J^ + O^- 



(71) 



The terms c^ 2 '™-' for n = 1,2,3,4 are determined exactly by NNLL soft-gluon resummation for 
arbitrary m t [4]. On the other hand, only parts of the c^ 2,0 - 1 coefficient are determined by the 
NNLL calculation, namely its //-dependence and the contribution from the product of NLO 
corrections to the hard and soft functions. To determine this coefficient exactly would require 
the massive hard and soft functions at NNLO in fixed-order perturbation theory, which are 
parts of the N 3 LL calculation. 

We now discuss to what extent the results from this paper can improve the NNLO approx- 
imation described above. For the pieces proportional to the Laplace-space logarithms there is 
no possible improvement to the exact results as a function of m t . However, it is a non-trivial 
check on the factorization formula (120]) that its expansion in fixed-order reproduces these re- 
sults in the limit m t /M — > 0. We have confirmed explicitly that this is the case, also for the 
rih pieces after converting the results for the massless hard and soft functions to a theory with 
five active flavors using (1641) . For the c^ 2 ' -* term, on the other hand, we can take the further 
step of determining exactly the terms enhanced by logarithms of M/m t . Moreover, since the 
NNLO fragmentation function is known, and all diagrammatic calculations needed to extract 
the contribution of the two-loop hard function to this term are in place, the only piece needed 
to fully determine this coefficient in the limit m t /M — >■ is the NNLO massless soft function. 
This is a much simpler calculation than that for generic top-quark mass, and once completed 
it will provide valuable insight into the uncertainties in approximate NNLO calculations based 
on NNLL resummation alone. 

Given results for the Laplace-space coefficients, we can obtain a result in momentum space 
by making the replacements in (1451) . The momentum-space coefficient then takes the form 



C (2) (z,M,m t ,cos9,n) = D 3 



ln d (l - z) 
1 - z 



+ D 7 



r ln 2 (l - z) 

1 - z 



ln(l - z) 



D r 



1 - z 



+ C 5(l-z)+R(z) 



(72) 



The coefficients Dq, . . . , D3 and Cq are functions of the variables M , m t , t\, and /i. The plus 
distributions Di are determined by the approximate NNLO formula in Laplace space; explicit 
results were given in [4]. These plus-distribution coefficients are exact, valid for generic values 
of the top-quark mass. To determine the delta-function coefficient exactly would also require 
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the unknown pieces of the c^ 2 ^ coefficient, in other words the NNLO hard and soft functions. 
Without these pieces, there is an ambiguity as to what to include in this term. Parts of these 
are directly related to those in Laplace transformed coefficient c^ 2 '°\ and parts are related to 
whether to include the results from inverting the Laplace transform via the replacements (145]) ; 
we will always specify our means of dealing with these ambiguities when giving numerical 
results later on. 

The decomposition fT72|) is the natural one for studying the NNLO correction within the 
framework of soft-gluon resummation. However, it is worth mentioning that another way to 
estimate higher-order corrections would be to instead focus on the leading terms in the m t /M 
expansion, as determined by factorization formula (1141) for the small-mass limit. Given NNLO 
fragmentation functions and the NLO massless scattering kernels CV,- for generic values of z, 
it would then be possible to study to what extent the leading terms in the soft limit reproduce 
the singular logarithmic corrections in the limit m t /M — > 0. While such a calculation would be 
interesting and valuable as a means of studying power corrections to soft-gluon resummation, 
it is clearly beyond the scope of this work to pursue this idea further. 

5 Numerical studies 

In this section we perform short numerical studies of our results. Since a detailed analysis 
of soft-gluon resummation for the pair-invariant mass distribution was carried out with the 
counting m t ~ M in jUE|, our main motivation is to study to what extent these results can 
be improved through the additional layer of resummation in the small-mass limit m t -C M. 

When arguing that resummation is required in a certain limit, a typical first step is to 
check to what extent the corrections at a given order in perturbation theory are related 
to the logarithmic pieces. If the logarithms account for the bulk of the corrections, it is 
an obvious improvement to use the resummed formula to include subsets of higher-order 
corrections related to them. Comparisons of NLO corrections in the z — > 1 limit with exact 
fixed-order results were performed in |U[5j|25]. It was shown there that the logarithmic plus 
distributions, determined by an NLL calculation, account for a bit more than half of the NLO 
corrections, while these logarithms plus the delta-function term, determined by an NNLL 
calculation, account for essentially all of it. Moreover, it was observed that the perturbative 
corrections at the scale fif = M are rather large at high values of the invariant mass. The 
main question we seek to answer in this section is whether this is due to the small-mass 
logarithms. If so, it would be necessary to supplement the phenomenological results of [5] 
with the small-mass resummation derived here. 

We address this question at NLO by isolating the terms which can give rise to large 
logarithms in m t /M. This is easily done by expanding the resummed formula (I6TJ1) to NLO 
in a s (fif), for the choice \ih = Us = M, fidh = Hds = Ht = iTit- The NLO corrections 
proportional to mass logarithms are determined by the NLL calculation, and when expressed 
in a theory with five active massive flavors, they read (normalized to the Born-level Laplace- 
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space coefficient c^] 




1 cusp,0 in 2 cusp,0 

m t 




9 \ A/f 2 / \ A/f 2 \ 1 

3 6 q «N h + 7 *j In ^ + (-A) + 7o ) In J 5(1 - | . (73) 

We have checked numerically that for /iy ~ M these small- mass logarithms make up only a 
small part of the NLO corrections in the z — > 1 limit (which in turn account for most of the 
full correction), even for values of M as high as 3—5 TeV at the LHC with y/s = 7 TeV. The 
NLL corrections determined by soft-gluon resummation for arbitrary m t , which include also 
plus distributions containing no small-mass logarithms, make up a much larger part of the 
exact NLO correction. We will therefore view small-mass resummation as supplementary to 
that, a means of improvement rather than a substitute. The point is that the mass logarithms 
generated at NLL in the small-mass limit are a subset of the NNLL corrections in soft gluon 
resummation for generic m t , so including them may be advantageous. We now study this 
statement in more detail. 

To do so, we define different approximations to the NLO corrections at the level of the 
Laplace-space coefficient, which reads 

2 M 2 ( 1 \ 

cW (N, M, m u cos M/) = E c (1 ' n) (M, m u h, N ) In" _^ + O ( - J . (74) 

n=0 " f 



Momentum space results are obtained through the replacement rules (I45I) . The coefficients 
c^ 1 ' 2 ) and c*- 1,1 ^ thus determine the plus-distribution coefficients. We will take the exact ex- 
pression for these, valid for arbitrary m t . We then distinguish between three approximations, 
which differ only in their treatment of c^ 1,0 ^: 

1. use no information, i.e. c^ 1 ' - 1 = 0; 

2. use the information from the NLO fragmentation function, the heavy-flavor matching 
coefficients, and a s decoupling in the massless hard and soft functions, thereby including 
all terms enhanced by In mt/M for \if ~ M; 

3. use the information from approximation 2 plus the constant pieces of the massless NLO 
hard and soft functions. 

These three approximations add progressively more information to the Laplace-space coeffi- 
cients: the first is NLL in soft gluon resummation, the second contains parts of the NNLL 
correction enhanced by small-mass logarithms, and the third includes the entire NNLL correc- 
tion but expanded in the limit mt/M — > 0. We show the NLO correction to the invariant mass 
distribution obtained from these three approximations in Table [TJ We also show the "exact" 
result obtained from the leading terms in the z — > 1 limit, for arbitrary top-quark mass. We 
see that the first approximation accounts for about 50-60% of the exact answer. The second 
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M = 500 GeV 


71 /r -i v r\r\ ^ T 

M = 1500 GeV 


M = 3000 GeV 


Approx. 1 


0.074 


1.04 x 10~ 4 


0.70 x 10~ 7 


Approx. 2 


0.085 


1.35 x 10" 4 


0.94 x 10~ 7 


Approx. 3 


0.126 


1.79 x 10" 4 


1.19 x 10~ 7 


Exact (z — > 1) 


0.154 


1.86 x 10" 4 


1.20 x 10~ 7 



Table 1: NLO corrections to the differential cross section da/dM (in pb/GeV) computed using 
m t = 172.5 GeV, /i/ = M and MSTW2008NNLO PDFs, at LHC with ^s = 7 TeV. 





M = 500 GeV 


M = 1500 GeV 


M = 3000 GeV 


Approx. A 


5.67 x 10~ 2 


1.22 x 10~ 4 


1.11 x 10~ 7 


Approx. B 


6.35 x 10~ 2 


1.40 x 10~ 4 


1.26 x 10~ 7 



Table 2: NNLO corrections to the differential cross section da/dM (in pb/GeV) computed 
using m t = 172.5 GeV, /i f = M and MSTW2008NNLO PDFs, at LHC with ^ = 7 TeV. 

approximation is a small improvement, and can account for up to 75% of the exact answer 
for high invariant mass. The third is a big improvement over the first two, especially at lower 
values of the invariant mass, and accounts for nearly the entire correction in the z — > 1 limit 
already at M=1500 GeV. 

We now discuss approximate NNLO corrections, and how they change upon including extra 
information from the small-mass limit. Following our NLO analysis, we take the Laplace- 
transformed coefficient (171 p as the fundamental object. As mentioned in the previous section, 
the terms c^ 2,n ^ for n = 1,2,3,4 are known exactly for arbitrary m f , since they are part 
of the NNLL calculation. We keep the full mass dependence of these terms in our NNLO 
approximation, and convert them to momentum-space results as in fj4"5j) . We then consider 
two approximations for the constant term c'- 2,0 -': 

A. use no information, i.e. c 1 - 2,0 ^ = 0; 

B. use the information from the NNLO fragmentation function, plus the rih terms aris- 
ing from a s -decoupling and the heavy-flavor coefficients, thereby including all terms 
enhanced by (up to two) powers of In mt/M for /i/ ~ M . 

The first of these is a pure NNLL calculation in soft gluon resummation for arbitrary m t , while 
the second is NNLL plus the part of the N 3 LL correction enhanced by small-mass logarithms 
(as well as some unenhanced terms associated with the fragmentation function, heavy-flavor 
coefficients, and a s -decoupling to second order). The results for the NNLO correction cor- 
responding to these choices are shown in Table |5J The two approximations are numerically 
rather close to one another. This shows that logarithms enhanced by powers of In m t /M do 
not lead to large corrections. Moreover, at higher values of the invariant mass the NNLO 
corrections can be even larger than the NLO ones. This motivates all-orders soft-gluon resum- 
mation instead of NNLO expansions, which was indeed the approach taken in [5]. However, 
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the only way to know the size of the missing N 3 LL corrections is to calculate them, and it will 
be interesting to return to this issue once they are known in the mt/M — > limit. 

6 Conclusions 

The pair invariant-mass distribution is an important observable for top-quark physics at 
hadron colliders. In this paper we set up a framework for dealing with potentially large 
perturbative corrections at high values of invariant mass. In particular, we gave explicit fac- 
torization and resummation formulas appropriate in the double soft (z — > 1) and small-mass 
(m t <C M) limit of the differential cross section, along with the ingredients needed to evalu- 
ate them to NNLL order. While many of the ideas and perturbative calculations needed to 
accomplish such a resummation were already available in the literature, this is the first time 
they have all been combined for a description of the invariant-mass distribution in top-quark 
pair production at hadron colliders. With small modifications the methods can also be used 
for the double soft and small-mass limit of single-particle inclusive observables such as the px 
or rapidity distribution of the top quark. 

We deferred a detailed phenomenological study of our results to future work. However, a 
short numerical study revealed the following features. First, even for invariant mass values as 
high as 3 TeV, it is not obvious that small-mass logarithms of the ratio m t /M are so large 
that they need to be resummed. On the other hand, already for values of the invariant mass 
of around 1.5 TeV, the leading terms in the m t /M — > limit, including the constant pieces 
as well as the logarithms, provide an excellent approximation to the full correction within the 
soft limit. We thus envision the main utility of the results obtained here as a means of adding 
parts of the N 3 LL corrections to soft-gluon resummation for arbitrary values of m t , as an 
expansion in m t /M. The missing piece of this analysis is the NNLO soft function for massless 
two-to-two scattering. Once completed, the calculation of this function will provide valuable 
information into the importance of higher-order corrections to the NLO+NNLL results for the 
invariant mass distribution presented in [5]. 
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A Small-mass limits of massive hard and soft functions 



Our derivation of the factorization formula (I18j) used as a starting point the small-mass fac- 
torization formula (]14p for the differential partonic cross section. In this appendix we briefly 
discuss the alternative derivation starting from the factorization formula (fTUj) and studying the 
properties of the massive hard and soft functions in the small-mass limit. For simplicity, we 
neglect contributions of heavy-quark loops, which give rise to additional terms proportional 
to rih taken into account by the heavy-flavor matching coefficient (TT91 

We begin by recalling that the massive hard matrix is related to UV-renormalized virtual 
corrections to color-decomposed amplitudes for two-to-two scattering. These IR-divergent 
quantities are rendered finite through multiplication by a renormalization matrix Z m . We 
define the bare and renormalized quantities according to 

limZ- 1 (e,M,m t ,t 1 , / u)|M(e,M,m 4 ,t 1 )) = \M reil (M, m t ,t u »)) , (75) 

where the quantity on the right is finite in the limit e — > 0. We write the analogous relation 
for massless amplitudes as 

\\mZ- x {e,MM^) \M{e,M,t x )) = \M ien (M, t u /i)> . (76) 

We can use these together with the relation between massive and massless scattering ampli- 
tudes in the small-mass limit [27] to derive a relation between the massless and massive hard 
functions. To do so, we first use 

\M{e,M,m t M)) = Z [q] (e, m t , \M(e,M,h)) (77) 

where Z\ q ] is the object given to NNLO in Eqs. (37) and (38) of [27]. Here and in the remainder 
of the section we have neglected terms which vanish in the limit m t —> 0. We then multiply 
both sides of (1771) by the massive renormalization factor, use (1751) to deduce that both sides 
are finite, and then observe that (17611 implies that 



Z[ q ](e,m t ,ti)Z m 1 (e,M,m t ,t 1 ,ti) = f(m u /j,)Z l (e, M, t u //) . (78) 

The function / is a scalar matching correction which is finite in the limit e — > 0. From this 
relation and the definition of the hard matrix in terms of the IR- renormalized color decomposed 
amplitudes, it then follows that 

H${M, m t , h, fj) = f 2 (m t , ^H^M, t uf i). (79) 

Since when neglecting terms proportional to nh the only m t dependence in the factorization 
formula for the partonic cross section is through the fragmentation function, it must also be 
true that 

S%(V§(l-z),m t ,t 1 , f j,)= (80) 
rx C D (m t ,ii)S D (m t (l - z),fi) C D (m t ,ii)S D {m t (l - z),n) 

- Z ), h ,„) . 
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Finally, we recall that the soft function is related to real gluon emission in the soft limit, a 
statement independent of whether the final-state quarks are massive or massless. Since all 
real emission contributions in the fragmentation function are associated with the function Sd, 
we expect 

f(m t , fx) = C D (m t , //) , (81) 

which would imply a simple relation between the massive and massless soft functions in (180|) 
through a double convolution with partonic shape functions. We have checked that (1801) and 
( IHTj) are satisfied for the NLO functions, and for the part of the NNLO functions determined 
by approximate NNLO formulas. In addition, we can use ( 18 ip as a consistency check between 
various NNLO calculations available in the literature, namely those for the renormalization 
factors in (I78p . that for the fragmentation function, and that for the shape function. We find 
nearly total agreement, with the exception of a piece related to the CaCf color factor. In 
particular, we find that direct evaluation of (T78|) to NNLO yields 

f(m t ,fi) = C D (m t ,fi)-4n 2 C A C F {^f . (82) 
Unfortunately, we have not been able to resolve the source of this discrepancy. 



B Convolutions of plus distributions 

Let us define the function 



3) 



2r(2?7) (1 - z) 1 -^ ' 

Representations of plus distributions can be obtained by taking derivatives of / with respect 
to rj and by subsequently expanding the result in the limit 77—7-0. One finds for example, with 
a test function g, 





f dz 


J 









f 


dz 


Jo 




dz 


"ln(l 




1 



dz5{l — z)g(z) 



dzf(z,rj)g(z] 



?7— 5-0 



1 - Z 



g(z) = d v I dzf(z,rj)g(z) 
1 



»?-s-0 



9(z) = -d n I dzf(z,r))g(z) 



7T 



r?— >0 



+ - / dzf(z,ri)g(z) 



54) 



T7— 5-0 



These representations are useful when one wants to calculate the convolution of plus dis- 
tributions. Consider for example the following convolution 



1 



1 - z 



1 



1 -z 



1 dx 

x 



1 — X 



1 — z/x 



55) 
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By employing the differential representation found above one finds 



1 -z 

The integration gives 
' dx 



1 - z 



d d 



dx 



x 



(86) 



??l->-0,7?2->'0 



X 



f(x,Vi)f (~'^J 



-2(j7i+r/ 2 )7i 



- 2 F X (2 Vu 2^, 2(771 + 772), 1 - *) (1 - z)- 1+2 ^ + ^ . 



57) 



Since we are only interested in the terms which are singular in the z — > 1 limit, for the purposes 
of our calculation one can actually set 2 -Pi — > 1- The limits for 7?j — > are then considerably 
easier to evaluate. To do so, one makes the following analytic continuation in the integral in 



(1 _ z yi+i(vi+m) _> (1 _ z) ~i+2{r n+m) + 8(1 z) 

2(771 + 772) 

One can then safely take the derivatives with respect 771 and 772 in Eq. f[8"B"j) and take the limit 
for vanishing 77, to obtain 



1-z 



1 - z 



ln(l -z) 
1 - z 



C 2 5(l - z) . 



(89) 



The same kind of procedure can be employed to calculate other convolutions. For example 



ln(l - z) 
1 - z 



1 - z 



1 12 



dl \ fUx 

2 



1 +C2) / —f(x,m)f[-,V2 



x 



ln'(l - z) 
1 - z 



1-z 



+ Cs5(l-z) 



(90) 



m(l - z) 
1-z 



ln(l - z) 
1 - z 



\n 6 (l - z) 
1 - z 



+ 



m->-o 



-2C2 



Mi - g) 

1 - z 



2Ca 



1 - z 



r?2->-0 

C(4) 



(91) 



C Anomalous dimensions 



The anomalous dimension 7 cusp , introduced in (I28f29p . has the following expansion in powers 
of a,: 



n. 



(92) 
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Completely analogous expansions hold for 7 9 , 7 9 , 7 s , and 7^. The coefficients of the expansion 
in (Ml are EH 



cusp 

7o 



cusp 

7i 



cusp 
72 



268 4vr 2 



9 



C 2 



^ 80 

C A - — T F n f 



490 536tt 2 44vr 4 88 _ 



f672 160vr 2 224 , 

1 C3 

27 27 3 s 



C F T F n 



I 



220 



64£ 



64 _o n 

27 T ^<- 



(93) 



The coefficients in the expansion of 7 9 and 7 9 up to 0(a 2 ) are [55|l56 



it = -3C f 



ll = C F \-- + 2tt 2 - 24C 3 + C F C A 



961 llTT 2 



54 6 



26C 3 +C F T F nj 



130 2tt 2 

IrT + IT 



(94) 



and USES 



7o = — g-CU+ 3^>n/ 

692 
~27~ ' 



lt = C. 



(95) 



The coefficients in the expansion of 7 5 up to 0(a 2 s ) are [TT } IT2" } 
7 5 = -2^f , 



7f = C f 



110 7T 5 



I8C3 cu 



8 2 

h -' 

27 9 



7T J T F Uf 



27 18 

The coefficients in the expansion of PDF anomalous dimensions up to 0(a 2 s ) are 
it = 3CV , 



(96) 



it 1 = c\ 



■ 2tt 2 +24C 3 +C F C. 



17 22tt 2 
"6 _ + _ 9 _ 



12( 3 \-C F T F n f I 



2 8tt 2 



3 9 



(97) 



and 



7o = y^A - -Tpra, 
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( y + 12Cs) - f C A T F n f - 4C F T F n f . 



for the gluon and quark PDFs respectively. 

Finally, we define expansion coefficients for the QCD (3 function as 



-2a K 



s 

in 



where to three-loop order we have 
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yC'l ~ ^C A T F n f - 4C F T F n f 
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(98) 
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